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1. INTRODUCTION 

The essential preference of the MESFET device is based on the high mobility of electrons in the 
channel compared to the MOSFET device. modelling transistors I-V characteristics is an essential of any 
circuit design [1], when can be replaced in many applications, the most important applications today are cell 
phone technology, optical fiber communications systems, high performance electronic test equipment, and 
military transmission applications [2], [3]. MESFET-GaAs (metal semiconductor field effect transistor) made 
from gallium arsinide which is the basis of fast digital integrated circuits and microwave circuits and 
indispensable for applications operating at frequencies >50GHz. Important efforts have been deployed to 
increase the power and frequency performance of MESFET-GaAs and MESFET-GaN [4], [5]. The structures 
for these MESFET’s are manufactured by two methods: ion-implantation and epitaxial grow [6], [7], the 
output characteristics of ion-implanted transistors are highly dependent on the quality of the substrate [6]-[8], 
for MOSFET, N20 is the thin insulator deposition procedure used in the dielectric polysilicon gate of thin 
film transistors to achieve SiO2 layers at temperatures below 150°C [9]. 

Several analytical models have been proposed, among his known model have found the model of 
Larson [10], the model of Curtice and Curtice-Ettenburg model [11], [12], the disadvantage of these models 
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is that they cannot simulate the dependence of the parameter V; on Vps and Vas, so they are not suitable for 
short-channel devices [12], [13], one of the problems of the Tajima et al. [14] is that it does not take into 
account the simulation of Gm and Ga, for the simulation of MESFET-DC characteristics, Rodriguez 
introduced an expression based on an expansion of the Curtice model depending on the parameter V+ on the 
Ips expression [13], the adjustment of the precision of the Curtice model at the level of the linear and 
saturation regions is improved by the Chalmers model [15], let’s look at the Dobes model [16] that we notice 
for n=2, we return to the Rodriguez model [17], the Materka-Kacprzak model [18] Seemed to be a bit more 
precise in the linear and saturation regions, but its accuracy degrades significantly when the device size is 
reduced [19], the Memon model simulates the characteristics of a MESFET-GaAs having a finite density of 
states at the Schottky barrier [20], by modifying the Ahmed model [17]. It has been shown that the Memon 
model can simulate the characteristics of the device[20], we therefore propose another new model which 
takes into account the interface states at the level of the Schottky barrier, which takes into account the 
simplicity of determining the different characteristics of the MESFET-GaAs device, such as the output 
transconductance, the output conductance of the MESFET device. 

To improve device performances and integration, device dimensions are drastically reduced. This 
results in several short channel effects that affect the MESFET behavior and functionning. Reducing the size 
of the device conducts to a sometimes too complicated modeling of the drain current of MESFET's devices. 
Several models of current expressions have been presented and developed to explain the behavior of the drain 
current in the linear and saturation regions of the I-V characteristics, a typical explanation for the MOSFETs 
of Benfdila et al. [21], [22] and this is found to be the same for MESFET transistors, these effects limit the 
device performances by reducing the conductance Gm and increasing transconductance Ga when the gate 
length gets smaller. In this work, we describe a GaAs-MESFET modelling of the MESFET-GaAs structure, 
later proposed an analytical model that approaches the simulation and then we compare the results obtained 
with the proposed model and thus with the Curtice and Rodriguez models [11 ]-[13], the simulation includes 
fundamental equations which take into account the potential and the density of charge carriers, the poisson’s 
equation, the transport equations, and the continuity equations [22], [23]. The impurity conditions of the 
trapped recombination and the various parameters that influence the static output characteristics of the 
MESFET device, the simulation includes the effect of temperature, to see the behavior of the device under 
the effect of the temperature settings on the drain current curves of the MESFET-GaAs device. 


2. DEVICE STRUCTURE AND TECHNOLOGY 

The basic structure of the MESFET-GaAs device taken in this simulation is shown in Figure 1. The 
deep zone is considered to be semi-insulating GaAs doped with a concentration of Na, the channel is doped is 
an N-type with Na concentration. The implanted MESFET consists of a minimal concentration semi- 
insulating GaAs substrate with an implanted channel of uniform profile and maximum concentration 
Nu=l.e!%cm °, the thickness of the channel is d=0.12um. The n* zones implanted under the source and drain 
electrodes with a thickness of 0.054m and they have a maximum concentration Na*=1.e!8cm>°. The gate 
length L=0.3um, the aluminum metal was selected for the Schottky contact, the working function of the 
aluminum is 4.28ev, the length of the source-gate and the length of the gate-drain Lg=Leq=0.25ym, the 
values of Tn and Tp have the value of 10 [24], [25], the saturation velocity v.=2.3xe°m/s [25], The substrate 
has a thickness of 0.3um, the total thickness of the device is 0.42um and its maximum length is 1.2um. The 
simulations describe the drain current, transconductance variations, the output-conductance at different bias 
voltage and the influence of the temperatures on drain current characteristics. These are key points discussed 
in this work regarding the GaAs-MESFET performances. 


Lgs Lg Lds 
mm + —_ 
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Figure 1. GaAs-MESFET device structure and geometry (after [26]) 


As well as the temperature is one of the essential parameters to take into account. Indeed, the 
temperature makes it possible to modify the performances of the device [23], [27], [28]. This work has 
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developed the goal of finding a new analytical model based on our numerical simulation and compared to the 
two different I-V nonlinear Curtice and Rodriguezmodels for the GaAs short-channel MESFET. The Figure! 
shows a two-dimensional of the GaAs-MESFET structure simulated in this work, the channel has a thickness 
of 0.12um and is doped Na=1.e'*cm%, the substrate has a thickness of 0.3um where the shallow residual 
impurities N.=1.e°cm > are compensated by Ny=1.e!°cm > traps. 

We have chosen the device parameters to fit the device design requirements for short channel 
MESFETs. These physical and geometrical parameters are typical for a short channel GaAs-MESFET device 
specifications and typical parameters of submicronic MESFET used in high speed digital and analog 
electronics are used in the proposed model are listed in the Table 1. 


Table1.Parameters model of the MESFET after ([29] and [30]) 








Simulation parameters Values Unites 
Nz Channel doping concentration 1.e'® Vem? 
Nat Source / drain doping concentration 1.e!8 I/cm? 
N, Substrate doping concentration, l.e? Iceni 
d Channel Thickness 0.12 um 
L; Metallurgical Gate length 0.3 um 
Lgs Gate-Source spacing length 0.2 um 
La Gate-Drain spacing length 0.2 um 
Voi Built-in Potential 0.87 V 
V,a Saturated value of electron drift velocity 2.3e5 cm/s 
Nr Substrate compensation trap l.e” Vem? 
Us Electron affinity of GaAs 4.07 eV 
Øm Gate metal work-function 4.28 eV 





2.1. Numeric simulation 

The modeling of the two-dimensional geometry of the grid short-channel MESFET including 
Classical semiconductor equations such as the Poisson’s (1) and electric potential defined in terms of the 
electric field (2), the simulation involved the solution of the Poisson’s, continuity equations, two nonlinear 
partial differential equations[3 1]. 


Pixy) 
A Woy) =-= (1) 


EGaAs 


Egy) = —VWey) (2) 


The expression of the current density in relation to electric field is given by: 


> > 


Jay = Oxy) Ex,y (3) 


The MESFET was considered of type-N, where the density of current I was not aken into account during 


= 


the whole simulation, which is evident for the total current density ise = jee ie + mA where A n dyn 
are the density of the electron currents in the x and y-directions respectively. 
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Joy = WrnEx + qD Vn (4) 


=>. 


Jya = qunnEy + qD, Vn (5) 


The electron concentration density can be written in the form n = nno exp( — qVbi/KT), where V»i defines 
the built-in voltage whose value is given in Table 1 and Dy, is the diffusion of the electrons constant is given 
by the Einstein relationship (6), and the current density total from the channel (7). 


Dn = hip (6) 
q 

> 2 oR dn dn 

Jtot = 9 (unn (E, F Ey) T (= + = Dn) (7) 
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The physical quantities involved in these equations are as follows: 

- Egyyis the electric field on x and y-directions 

- ‘Poy is the total variation of the potential in the space charge zone 

-  Jx is the conduction current density in x-direction 

- Oy), Po.yyare the conductivity and the charge density in the channel respectively 

Solving the Poisson equation to extract the expressions the electrostatic potential wy. 

Determination from (2) the components E, and Ey of the electric field E &,y) in the channel. Integration then 
determines the current J, at the source and drain electrodes, in order to establish the analytical expression of 
the current Ia. The current Ia at the terminals of the device is generally calculated by integrating the total 
current density Jt over an appropriate surface in each contact zone (8). 


Ig = [Jenas (8) 


The simulation contains the implementation of the Poisson’s, the current continuity equations whose 
flowchart is presented in appendix-A .1. The Dirichlet (fixed values) and Neumann (zero derivative values) 
boundary conditions for the potential are applied to the contacts and to the surface, respectively for the device. 


2.2. Classical current, conductance and transconductance models 

The classical models are investigated during decades and have been improved to match submicron 
MESFETs as done earlier for submicron MOSFETs. This is targeting their use in high speed digital circuits 
and an eventual substitution or completion of advanced CMOS technology. However, the temperature effects 
were not deeply studied. In this work we will introduce these effects and hence introduce a new model. The 
basic current equations of drain current expressions of different nonlinear MESFET of the Curtice and 
Rodriguez model respectively are given by the two expressions (9) and (10) respectively: 


Ins = Bes — Vz)? tanh(aVp5)(1 + AVps) (9) 
Ins = B(Ves- Vi — YVps)? tanh(aVps)(1 + AVps) (10) 


Where £ is the transconductance, a simulate the linear region dependence, À simulates the dependency of Ips 
in the saturation zone on Vps, y simulates the threshold voltage dependency on the Vps, V: the threshold 
voltage, the V:and AV; are defined by [20], [32]: 


In qNad* 


V, 
i 2E; 





D, (11) 


Where q is electron charge, d is the channel thickness, €s is permittivity of GaAs, N4 the channel doping and 
@, is Schottky barrier height. 
4d 


AV.= > V, 
a (12) 


AV; is the threshold voltage due to the submicron structure and Ly is the gate length [33], [34]. 
2.2. Conductance and transconductance equations 


The MESFET transconductance is defined by a ratio of variation of the current Ips by the voltage 
Vas (13) and also the conductance obtained by the variation of Ips by the voltage Vps (14) [20]. 





Gq = Vg Ves = Cst (13) 
Gm = Ves Vos = Cst (14) 


The expression of the conductance Gaand of the transconductance Gm, respectively of the Curtice model 
obtained by directing (9) according to Vps and Vas respectively are (15) and (16) respectively. 


i 2a F: À 
PS [sinh(2aVps) 1 + YVps 





Ga = (15) 
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Chas ee (16) 


Vps-Vt 


Concerning the expressions of the conductance and the transconductance of the Rodriguez model are 
obtained by directing the (10) following Vps and Vas respectively are (17) and (18): 


G, =I ee a a a AN (17) 
a DS [Vos —Vi—YWos sinh (@Vps) 1+ YVps 
2I 
Gm = a eaa (18) 


Vos — Vte — YVps 


3. PROPOSED MODEL 

Several modeling approaches based on physical characteristics have been used to perfect and 
improve the electrical and thermal performance of GaAs MESFET transistors. often have a finite state 
interface density and the ideal device pinch. Certain physical parameters such as the barrier height, the 
density of the interface states, generate problems of the interface states, substrate compensation traps and 
non-homogeneous Schottky barrier, the model includes the voltage V, and the Schottky potential contact at 
the same time given in (19). In submicron GaAs-MESFETs, in order to study the submicron MESFET, we 
have simulated and drawn the device current-voltage (I-V) characteristics for different drain and gate 
voltages for temperatures ranging from 200K to 600K steped by 100K. 

Assuming Na=1e!6 1/cm? and d=0.12um and knowing the value of V, using (11), which depends on 
the Schottky potential ®,, the potential ®, is determined by the difference of the work-function of the metal 
Ọm and the electron affinity 7, of the GaAs semiconductor which is equal to ®,,-7,, the AV; is obtained by 
(20), which represents the shift in the threshold voltage due to the submicron structure of the MESFET whose 
gate length L,=0.3um for the devices, whose parameters are and represented in Table 1. The table shows that 
the device having L;=0.3um it is with this amplitude that the transistor is classified as submicron, increasing 
the channel thickness d, the performance of the component deteriorates. The basic current equation used in 
this study are given by the following current-voltage expression: 


T 
Ins = Bae (A)? tanh(aVps)(1 + AVps) (19) 
0 


Where £ simulate the transconductance, à simulates the dependence of Ips on Vps in the saturation zone, V, 
the threshold voltage, a simulate the dependency of a linear zone, y simulates the dependence of threshold 
voltage on Vps, To represents the ambient temperature, T the temperature desired. 

The expressions of the output conductance Gy and the output transconductance Gm in small-signal of 
the proposed model are deduced therefrom the expression of the current Ips given in (19), with regard to Vps 
and Vas respectively are given by: 


T 
Ga = 3B pA lB — C]DE + aBFAPE + ABARE (20) 
0 
G 3B g A? B DE (21) 
m = 3B — — 
To Y 


The indices A, B, C, D, E, F are used for the simplification of the proposed model whose associated 
quantities for each index are given by: 





A= (1 —AVp5)B (22) 
B = J (Ves — V; — AV; — Vos) (23) 
C=y oe (24) 
D = tanh(a@Vps) (25) 
E = 1 +AVps (26) 
F = 1 — tank? (a Vps) (27) 
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4. RESULTS AND DISCUSSIONS 
4.1. Device DC characteristics 

The simulations are performed on the GaAs-MESFET considered in the structure of Figure | and 
the curves obtained underline the performance of devices. The I-V characteristics ID=f (Vps, Vas) for a 


device containing variations of the different physical and geometrical parameters of the proposed structure 
are presented in Figure 2. 
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Figure 2. Characteristics of the GaAs-MESFET at room temperature obtained by simulation, (a) Ips- Vps 
charactéristics with a variety of values Vas, the drain voltage from 0 to 2V by step 0,05V, 
(b) Ips-Vas charactéristics at different values of Vps 


The significant dependence of the device on the temperature shown in Figure 3 and examines the 
influence of the temperature parameter on the component's activity, is clearly visible on the transfer 
characteristics I-V=f(T) as it is indicated for 300K the current Ips=0.9mA for Vgs=0V, on the other hand for 
400K the current has increased about 0.05mA, the temperature changes the curved shapes in the linear region. 
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Figure 3. Characteristics Ips-Vps obtained by simulation with a variety of values of temperature, (a) Different 
values of bias Vas polarization, (b) Fixed Vas to 0V 


The Figures 4(a and b) show the characteristics Ips=f (Vps, Vas) for the different gate and drain 
voltage values, respectively, the results obtained by the simulation of the proposed 2D structure are compared 
with the proposed model, the results of the model are similar and reasonable by comparing them with 
Rodriguez, Curtice and the simulation, which brings us back to the validation of our model. 

The comparison shows a low average error between the simulation and the proposed model, 
however, the precision obtained is acceptable in the manufacture of electronic devices, the average errors of 
the different curves illustrated in Figures 4(a and b) are represented in the Table 2. 
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Figure 4. Comparison of features of GaAs-MESFET between the simulation and the model proposed, 


(a) Comparison of the Ips-Vps obtained by simulation and proposed model with a variety of values Vas, 
(b) Comparison of the Ips-Vcs obtained by simulation and proposed model with at different values of Vps 


Table 2. Comparison of the RMS errors and average error of GaAs-MESFETs 














RMS errors at different Vgs results of RMS errors at different Vps results 
Different Vas Figure 4(a) Different Vps of Figure 4(b) 
Simulation Proposed model Simulation Proposed model 
Vos=0V 0.13 0,16 Vps=0.2V 0.21 0.24 
Vos=-0.3V 0.24 0.27 Vps=0.4V 0.18 0.20 
Vos=-0.5V 0.37 0.41 Vps=0.6V 0.14 0.14 
Vos=-0.7V 0.61 0.64 Vps=0.8V 0.08 0.07 
Average Error 0.34 0.37 Average Error 0.15 0.16 





The RMS and average error values were calculated (28) for the different bias voltages Vos, Vps 
respectively for the models considered as indicated in Table 2 respectively, the precision of the model is 
evaluated by referring to the RMS error values as a function of both Vas and Vps, observing that the mean 
error of the model predicted with a very small deviation from the mean error of the results obtained by 
simulation in the two cases. To calculate the relative error where the error in percentage in this study of the 
different curves obtained during this work using: 





wins — Yi)? 


N * 100% (28) 


RMS(Error) = 


The RMSerror values as shown in Table 3 for Figures 5(a and b), the RMS error values for the 
various Vas values obtained by the simulation lower than that obtained by the Curtice and Rodriguez model 
respectively and also note the Average error values are also lower. In Figures 6(a and b), show that the 
suggested model is exhaustive and that it can simulate the output characteristics of the GaAs MESFET under 
different grid conditions, in comparison with already existing models. 


Table 3. Comparison of the RMS errors and average error of GaAs-MESFETs with a Curtice and Rodriguez 
models 
RMS errors for different Vos compared with simulation and proposed model 











es Rodriguez Model ___CurticeModel__ By Simulation Proposed Model 
Vos=0V 0.36 0.36 0.13 0.16 
Voes=-0.3V 0.37 0.38 0.24 0.27 
Vos=-0.5V 0.60 0.59 0.37 0.41 
Vos=-0.7V 0.60 0.62 0.61 0.64 
Average Error 0.48 0.49 0.34 0.37 
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Figure 5. Comparison of MESFET characteristics, (a) Ips-Vps comparison between the simulation and 
Curtice model, (b) Ips-Vps comparison between the simulation and Rodriguez model 


The RMS error values of Figures 6(a and b) are also presented in the same table, the results obtained 
by the suggested model are lower than those of the two models Curtice and Rodriguez respectively and also 
for the values Average error, including the results obtained by the proposed model, are lower than those of 
the two models, respectively Curtice and Rodriguez, as well as for the average error values. 
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Figure 6. Comparison the Ips-Vps characteristics of MESFET, (a) Comparison between the proposed model 
and Curtice model, (b) Comparison between the proposed model and Rodriguez model 


The parameters defined during the simulation of the proposed model are the following: a=1.2; 
B=0.665; y=—0.02; 2=0.001, in contrast to models Curtice and Rodriguez are: a=1.6008; f=0.345; 
y=— 0.006; A=0.17. 


4.2. Conductance and output transconductance 

The values of Ga and Gm for the device are presented in Figures 7(a and b), 8(a and b) for all Vps 
and Vas apply and different temperatures, comparing the simulation of the device and the proposed model, 
shows a good agreement that confirms the validity of the proposed model, we may assume that it is possible 
to simulate a larger range of MESFET's with greater precision with the proposed model. 

Figures 8(a and b), show the output conductance and transconductance of the drain, respectively as a 
function of Vps for the values applied to the gate and as a function of Vas for the voltages imposed on the 
drain, respectively in comparison between the proposed model and the simulation of the device, let's look 
model [17] does not take into account the effects of gate polarization Vps, on output conductance Gg, device 
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operation in the saturation zone were avoid in Islam et al. [35] made a modification to this model, showed 
that the effect of bias voltages (Vps-Vcs) on the output conductance of the built-in device. The proposed 
model shows an excellent precision concerning the output conductance presented in Figure 8(a), the 
transconductance represented in Figure 8(b) in comparison between the simulation and the proposed model 
show an almost similar, which brings us back to justify and validate this model to easily determine the 
different characteristics of the MESFET-GaAs device. 
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Figure 7.The transconductance, the conductance at ambient temperature obtained by simulation of MESFET- 


GaAs, (a) The output-transconductance Gm (Vos) at different Vps, (b) The output-conductance Ga(Vps) at 
different Vas 








T y T i T y T r T 
1000 H=- Benem jemma 
a =1.2; B =0:665 












—— Proposed model |.. —— Simulation 


Proposed Model 


1 
S 
œ 
Q 
© 





simulation 





i i Vas= 0V E 
y =—0.02; A = 0.001 Ss : Vps= 0.4V 
800 H-4 iene cet ate pense Ves= -0,3V + Vps= 0.6V 
Vas= -0.5V 600 Vpg= 0.8V 
Ves= -0.7V 








600 








400 leaa EA EE i 
400 H 


Output-Conductance Gd (uS) 





Output Trans-conductance Gm (uS) 
































o 1 2 3 4 5 20 
Drain Voltage (V) 


(a) (b) 


Figure 8. Characteristics of the conductance and the transconductance of the GaAs-MESFET device for the 
ambient temperature in comparison between the model suggested and the model of simulation, (a) The 
output-conductance Ga(Vps) at different gate voltage, (b) The output-transconductance Gm (Vas) at different 
drain voltage 


5. CONCLUSION 

All of this work contributes to the implementation of a simulation taking into account the influence 
of important parameters when there may be an influence on the characteristics of the MESFET GaAs 
transistor. Our study has been processed and written to determine the output IV, conductance and 
transconductance characteristics of the GaAs MESFET, to simulate the I-V features of the GaAs MESFET 
device, a completely new model has been developed for a short submicron channel. The results of the model 
established show a correlation well to our simulation results, the efficiency of this model compared with the 
models available in the literature, the RMS error and the average error were calculated and showing 
reasonable accuracy of the model, the model could be a useful tool for the design and characterization of 





Temperature dependent analytical model for submicron GaAs-MESFET (Mohamed Djouder) 


1280 O ISSN:2302-9285 


integrated circuits, it has also been shown that the new proposal is a valid model capable of obtaining all 
current-voltage GaAs-MESFET, conductance and transconductance characteristics, the effect of temperature 
on the output characteristics of the device was discussed. it has been shown that the temperature is one of the 
essential parameters to be taken into account. In fact, the temperature makes it possible to modify the 
performance of the components. As indicated, this type of device can operate at high temperatures appears in 
all the figures obtained as a function of the temperatures. The effect of these parameters on performance has 
been observed and presented, its results compare favorably with a submicron GaAs MESFET studied by 
several authors. 


APPENDIX -A 

A-1; A MATLAB script is used in the modeling of the device which takes into account the various geometric 
and physical parameters of which the main program flowchart is illustrated in Figure A-1. The simulation 
code is developed using the MATLAB software which consists in modeling and simulating the device by the 
finite element method for the prediction of the characteristics of the MESFET sub-micron transistor, the 
geometry of which we have taken in Figure 1 with the parameters shown in the Table 1. The simulation 
program can be used on the one hand for the design of field effect devices, the parametric configuration can 
be carried out to allow on the other hand, by calculating the various curves of the device, the simulation also 
takes take into account the effects of temperature variation as an important parameter that should not be 
neglected as long as the temperature changes and affects the correct operation of the device, the proposed 
model of the GaAs MESFET is developed to obtain a comparison between the simulation and the best-known 
Curtice and Rodriguez models, this model is implanted on the MATLAB, the model is presented in an 
exhaustive way where the current, transconductance and conductance have been presented. 
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Fig. A-1. Main flowchart of the computer program implemented under 
Matlab to simulate the differents characteristics 


A-2; The iterative process of optimizing the new analytical model proposed, including the optimal parameters 
corresponding to the results of our simulation and those of the Curtice and Rodriguez models were defined after 
the comparison of the data and the validation of the proposed model. The simulation steps model are: 1) 
definition of values: a, p, y, 4, a, Le; 2) calculate from V, AV, p; 3) writing of the current equation (Ips), 
conductance (Ga) and transconductance (Gm) equations; 4) test program; 5) comparison and validation. 
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